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Abstract 

The functioning and efficiency of natural photosynthetic complexes is strongly influenced by their 
embedding in a noisy protein environment, which can even serve to enhance the transport efficiency. 
Interactions with the environment induce fluctuations of the transition energies of and interactions 
between the chlorophyll molecules, and due to the fact that different fluctuations will partially 
be caused by the same environmental factors, correlations between the various fluctuations will 
occur. We argue that fluctuations of the interactions should in general not be neglected, as these 
have a considerable impact on population transfer rates, decoherence rates and the efficiency of 
photosynthetic complexes. Furthermore, while correlations between transition energy fluctuations 
have been studied, we provide the first quantitative study of the effect of correlations between 
interaction fluctuations and transition energy fluctuations, and of correlations between the various 
interaction fluctuations. It is shown that these additional correlations typically lead to changes in 
interchromophore transfer rates, population oscillations and can lead to a limited enhancement of 
the light harvesting efficiency. 
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I. INTRODUCTION 



The harvesting of sunUght by organisms such as plants and bacteria plays a crucial role in 
life on earth. Nature has optimized the photosynthesis process to a great extent, leading to 
a high efficiency in both the capturing of light and the subsequent energy transport within 
the organism's light harvesting complexes.-"- Pigment molecules in these light harvesting 
structures absorb the light, after which the excitation energy is eventually funneled towards 
the reaction center, where charge separation takes place and the energy can be utilized for 
biochemical purposes. A closer understanding of the functioning of these biological systems 
and what features are crucial to their high efficiency may provide further insight in how to 
optimize synthetic light harvesting structures. 

Recent experimental studies have shown that quantum coherence in light harvesting com- 
plexes persists on surprisingly long timescales.-"- While excitation energy transport (EET) 
in such molecular systems has been a subject of study for decades,—"— the aforementioned 
experimental observation of long-lived coherence despite its occurrence in noisy environ- 
ments at room temperature has spurred further theoretical research. Naively, one might 
think that coupling the electronic excitations to a dissipative environment will lead to rapid 
decoherence and a loss of excitation energy. However, it has been shown that environment- 
induced dephasing and noise may actually enhance the transport efficiency of excitations in 
such structures.—"— The rationale behind this lies in the fact that interaction with degrees 
of freedom in the bath can induce transitions of the excitation that are directed towards the 
reaction center. Recent studies have shed further light on the mechanisms that allow for 
optimal efficiency of the excitation, suggesting an intricate functioning where aspects such 
as the non-Markovian nature of the bath interactions^^i"— and correlations in the bath^^^^i"— 
may be exploited to achieve further optimization of the light harvesting complex. 

The interaction of the chromophores with their local environments leads to fluctuations 
of the transition energies and interactions. More speciflcally, changes in the environment 
will induce fluctuating local electrical fields. The resultant Stark shifts in the transition 
energies and transition dipoles of the chromophores will thus also fluctuate; in addition, 
environment-induced changes in relative position or orientation of the chromophores will 
also cause corresponding fluctuations in the interactions.— A number of studies have been 
performed where only the effect of transition energy fluctuations has been included, as 
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well as possible spatial and temporal correlations between these transition energy fluctua- 
tions.—!^"— This is commonly done in the framework of the extended Haken-Strobl-Reineker 
model,— which provides an attractive approach due to its simplicity and tractability. The 
original Haken-Strobl-Reineker model allows for a treatment of both transition energy and 
interaction fluctuations, and while it does not include correlations and inherently assumes a 
high temperature, the formalism can straightforwardly be generalized to remove the former 
drawback. Fluctuations in the interactions, in contrast to transition energy fluctuations, 
have typically not been accounted for in studies of photosynthetic complexes, despite the 
fact that these will generally occur as well. 

Fluctuations in transition energies and interactions are induced by variations in the local 
environments of the chromophores, and environmental variations around one chromophore 
will generate both fluctuations in its transition energy and in its interaction with the other 
chromophores. In addition, different chromophores may share part of the environment; both 
these arguments imply that the generated fluctuations will in general be correlated. While 
correlations between transition energy fluctuations have been studied,—*^"— there should 
generally also be correlations of transition energy fluctuations and interaction fluctuations, 
and between different interaction fluctuations. In previous studies, such correlations have 
commonly been discarded, even though there is no a priori reason to do so. A number 
of recent studies have also suggested the occurrence of such fluctuations in photosynthetic 
complexes.—"— Moreover, fluctuations of the various intermolecular interactions are usually 
not accounted for, although also these are typically induced by changes in the local en- 
vironment. We present a study of the effects of fluctuations of the interactions and their 
possible correlations with the various transition energy and other interaction fluctuations, 
and provide a numerical analysis of their measurable effects. When applied to natural pho- 
tosynthetic complexes, such as the Fenna-Matthews-Olson (FMO) complex,—"— it is shown 
that interaction fluctuations can have a considerable impact on the excitation dynamics and 
the efficiency, and these effects should generally not be neglected. Furthermore, we show 
that additional correlations can lead to enhanced oscillations of the exciton populations and 
modifications of the transfer and dephasing rates, and we quantify the dependence of the 
trapping efficiency on the correlations and initial conditions. 

The paper is structured as follows. In Sec. [Tll we provide the theoretical background of 
our study: in Sec. Ill At we introduce the Haken-Strobl-Reineker formalism that is used to 
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describe the effect of environmentally induced fluctuations and their possible correlations 
on the excitation dynamics, Sec. Ill Bl discusses the possible correlations, and in Sec. Ill CI 
we introduce the Fenna-Matthews-Olson photosynthetic complex that we will apply our 
theory to. In Sec. IIIIl we show the results of application of our theory to FMO: Sec. IIII Al 
elucidates the role of uncorrelated interaction fluctuations. Sec. IIII Bl focuses on the effect of 
correlations between fluctuations in transition energies and interactions, while subsequently 
in Sec. IIII CI we investigate the effect of correlated interaction fluctuations. Our conclusions 
are presented in Sec. [IVl 

II. THEORY 

A. The Haken-Strobl-Reineker model 

Upon absorption of the solar light and energy transfer into the FMO complex, an elec- 
tronic excited state will be created in the bacteriochlorophylls (BChl's). Due to the strong 
interchromophore interactions, this excited state (exciton) will be delocalized over a num- 
ber of BChl's. To describe the excitation and its dynamics, we employ the Frenkel exciton 
Hamiltonian,— 



h = Hs + hr + Hs-r = J2^^ i^) ^^i + '^^^ i'^) ^"^1 + a ^i^l^i + 

n n,m^n q 

where \n) describes a state where chromophore n is in its excited state while all others are 
in their ground states, Jnm is the interaction between chromophores n and m, and we have a 
bath of harmonic modes labeled by q. The first two terms are the system Hamiltonian Hs, 
the third term corresponds to the bath Hamiltonian Hr, and the final term describes the 
coupling between the system and the bath Hs-r, which we will keep unspecified for now. 
In our case, the system consists of the BChl's within the photosynthetic complex, while the 
environment includes any other degrees of freedom the excitations can couple to, such as 
the protein surroundings. To adequately describe the environmental effects, it is necessary 
to work in the density matrix formalism.—"— The density matrix corresponding to the wave 
function \ip{t)) is given by p{t) = \ip{t)) {ip{t)\, and evolves according to the Liouville-von 
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Neumann equation {h has been set to unity), 

t^ = [H{t)rp{t)]^C~p{t). (2) 

Here, we have defined the Liouville superoperator as £ = [H, ...]. The above time evolution 
equation of the density matrix is equivalent to the time evolution of the wave function 
\ip{t)) as given by the Schrodinger equation. At this point, we switch to the interaction 
picture, where the time evolution induced by the Hamiltonian terms Hs + Hfi, which are 
large compared to the system-bath interaction Hs-r, is explicitly removed from the time 
evolution of all operators A(t), 

Aj{t) = e'iHs+Hu)t^^^^^-i{Hs+Hn)t_ (3) 

We are primarily interested in the evolution of the system degrees of freedom, which we 
obtain by taking the trace over the bath degrees of freedom. This leaves us with the reduced 
density matrix p(t) in the interaction picture, which is the quantity we will consider from 
this point on when referring to the density matrix. 

The approach we take in treating interactions with the environment is based on a method 
first introduced by Haken, Strobl and Reineker.— In this approach, one models the bath- 
induced fluctuations as classical Gaussian Markov processes. While the original Haken- 
Strobl-Reineker (HSR) model assumed uncorrelated site energy fluctuations, one can extend 
the methodology to allow for correlations between the various fluctuations; we refer to Ref. 
41 for the details. The Hamiltonian can then be written as 



H = ^En \n) {n\ + ^ Jnm \n) (m| + ^ Vnmit) \n) {m\ . (4) 

n n,m,y^n n,m 

The terms Vnm{t) are the environmentally induced fluctuations of the various system Hamil- 
tonian matrix elements. Note that the Hermiticity of the Hamiltonian implies Vnm{t) = 
V^„{t). The averages {Vnmit)) can be set to zero without loss of generality, and since the 
fluctuations are assumed to be Gaussian Markov processes, the problem is fully defined when 
the correlation functions (Kim(^)Ki'm'(^')) known. The effect of the environment is now 
fully encoded in the bath correlation functions Cnmn'm'i'T') = (Kim(T)V^'m'(O)). 
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B. Correlation functions 



To proceed, we make the white noise assumption, where it is assumed that the bath relaxes 
on a time scale that is short compared to the exciton dynamics. The time dependence of 
the bath correlation functions is taken as a 5-function, 

Cnmn'm'iT) = (Km. (t) K'm' (0) ) = 'ynmn'm'S{T) . (5) 

Substitution into the time evolution equation of the density matrix, Eq. |2l yields 



dpnm{t) 

= —I 



{t) + ''Im'mnn' Pn' mm'n' Pn 



dt 

n'm' 

(6) 

At this stage, it is worthwhile to consider the correlation matrix 7 in more detail. First 
of all, the elements 'jnnnn and 7„mnm(^ 7^ are simply the variances of respectively the 
transition energy fluctuations and the interaction fluctuations. All other elements of 7 
correspond to correlations between different fluctuations. Secondly, the correlation matrix 
fulfills a number of symmetries, due to the requirement of a Hermitian Hamiltonian and a 
set of trivial index permutation invariances, 

7 abed = Icdah = 76acd = labdc- (7) 

We observe that there are generally three types of correlations, namely between fluctuations 
in 

• excitation energies (a = 6, c = d), 

• excitation energies and interactions {a = b,c ^ d or vice versa), 

• interactions {a ^ b,c d). 

A number of studies have included the flrst type of correlation, i.e. spatially correlated 
transition energy fluctuations, which have been shown to lead to appreciable changes in the 
transfer rates and efficiency of light-harvesting systems.— >2Si"— However, the second type of 
correlation has only been studied in the dimer,— while the third type of correlation has 
to the best of our knowledge not been included at all in studies of EET in photosynthetic 
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complexes. The central concern of this study is to provide an understanding of the relevance 
of such correlations to the time evolution of the various density matrix elements and the 
overall quantum efficiency. 

Besides the above symmetry considerations, additional restrictions apply to the various 
correlations in order to describe a physical system.— This becomes clear when one considers 
the expectation value of the difference between two fluctuation elements squared, which 
should obviously give a positive result, 

By considering such inequalities for various values of the indices, one finds a number of 
consistency conditions that need to be fulfilled, 

'fnmnm ^ 
2 I'fnmnm'l ^ 'fnmnm ~l~ ''inm'nm' ^ 

where the indices n,m and m' are allowed to be equal. In addition, the magnitude of the 
cross-correlation between two fiuctuating matrix elements is limited by the magnitude of 
the original fluctuations. It is straightforward to show that 

'^nmn'm' — y/'ynmnm'^n'm'n'm' ; 

where the equality holds when the fluctuations of the Hamiltonian matrix elements Vnm and 
Vn'm' are fully (anti-)correlated. These limitations on the allowed fluctuation correlations 
should be kept in mind when considering the effect of including the various additional 
correlations. 

C. The Fenna-Matthews-Olson complex and light harvesting efficiency 

The model light harvesting system for which we want to quantitatively probe the ef- 
fects of correlations in the environmentally induced fluctuations is the Fenna-Matthews- 
Olson (FMO) complex. FMO is the photosynthetic complex of the green sulfur bac- 
terium Chlorobium Tepidum, consisting of three weakly coupled, identical subunits of seven 
bacteriochlorophyll-a (BChl) molecules each.-"-i2ii22i2iii2iii^ Within such a subunit, after ab- 
sorption of the incoming light, the excitation energy is transferred to the reaction center 
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(11) 



where charge separation can take place. Due to their proximity to the primary hght har- 
vesting antennae, energy typically flows into the system at BChl's 1 and 6.— 

We model a subunit of the FMO complex by the following experimental Hamiltonian in 
the site basis (in units of cm~^),— 
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The corresponding exciton states, which we label in order of increasing energy, are given 
in Appendix lAl Of particular relevance for our numerical results are the exciton states s = 3 
and s = 7 with energies of respectively E3 = 224 cm~^ and Ej = 480 cm~^, which are the 
exciton states mostly associated with BChl's 1 and 2, and the lowest energy exciton state 
s = 1 with energy Ei = —24 cm~^ which is localized mostly on BChl 3, from where energy 
will be trapped to the reaction center. 

In order to describe the influence of the various correlations on the light harvesting 
efficiency, we introduce two competing decay channels.— 1^1^ Exciton decay, leading to an 
irreversible loss of the absorbed energy, is described by adding a decay term 

4t^^ = {ki'^ + ki^) /2, (13) 

where kli"^ is the exciton decay rate at chromophore n. In the FMO complex, energy is 
transferred to the reaction center, and we model capture of the excitation by the trap 
through a trapping term 

4r^^ = (4*^ + fc^^) /2, (14) 

where kn^ is the exciton trapping rate at chromophore n. In the FMO complex, the reaction 
center is located close to BChl 3, and we will only include trapping from that particular 
chromophore, kn^ = k^^">5n,i- 

The efficiency of the process can be defined as the branching ratio between energy trapped 
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at the reaction center and energy lost through exciton decay,- 



18^44.45 



Q = 




n 



(15) 




n 




n 



where r, 



n 



Jp°° dtpnnif) is the mean residence time at site n. 



III. NUMERICAL RESULTS FOR THE FMO COMPLEX 

As stated before, the FMO complex consists of seven coupled BChl's. In the follow- 
ing sections, we will first discuss the effects of uncorrelated interaction fluctuations, and 
subsequently proceed with a study of the effect of correlations between the various possi- 
ble fluctuations. Since many of the interactions are already small to begin with, we only 
consider fluctuations in the strongest interactions. Likewise, we only consider correlations 
involving the strongest interactions, i.e. interactions between BChl's of a magnitude exceed- 
ing 25 cm~^, which are shown in Fig. [H In addition, we only include correlations between 
site energy fluctuations of site n and interactions involving that same site, and similarly, 
correlations between interaction fluctuations that have one site in common. This choice is 
motivated by the fact that correlations between fluctuations are predominantly caused by 
shared local environments. 

The numerical calculations reported in the upcoming sections concern energy flowing into 
the FMO complex at BChl 1. As we label the exciton states s in order of increasing energy, 
this corresponds to initially exciting exciton states s = 3 and s = 7. One of our goals is to in- 
vestigate the effect of including energy-interaction fluctuation correlations on the experimen- 
tally observed oscillations of the exciton populations.- Therefore, we choose an initial excita- 
tion where exciton coherence is already initially present, = (|s = 3) + e'^ |s = 7)) / -\/2. 
Obviously, different values of the mixing angle 6 correspond to different initial populations 
Pn of BChl's 1 and 2. In particular, ^ = imphes Pi = 0.08 and Ps = 0.89, 9 = tt/2 yields 
p^^ p^ = 0.49, and for ^ = vr we have Pi = 0.90 and P2 = 0.08. 

Additional calculations have been performed for energy flowing in from the other side, 
that is, from BChl 6. This corresponds (mostly) to initially exciting exciton states s = 5 
and s = 6. The results do not change qualitatively, except for a considerable decrease in 
population oscillation frequency due to the smaller energy difference between exciton states 
s = 5 and s = 6 as compared to exciton states s = 3 and s = 7. For future reference, the 
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FIG. 1: A schematic view of the FMO complex, with the strongest inter-BChl interactions denoted 
by arrows. The excitation wih typically enter the complex at either BChl 1 or BChl 6, and will 
flow towards the reaction center which is most closely associated with BChl 3. 

exciton that most strongly overlaps with the trap state BChl 3 is exciton state s = 1. 

The parameters considered in our simulations have previously been used in Ref. Q to 
model the FMO system. We take all exciton decay rates equal at a value k^"^ = k^'^^ = 1 ns~^, 
while only trapping from BChl 3 occurs, with a trapping rate fc^*-* = 1 ps~^. Furthermore, 
we take 70 = 7nnnn. = 94 cm~^ for the transition energy fluctuations; note that there is a 
factor of 2 difference in the definition of '~fnnnn as compared to the quantity F in Ref. Q. 

A. Uncorrelated interaction fluctuations 

The effect of interaction fluctuations can be investigated straightforwardly with the cur- 
rent formalism, by introducing nonzero values for the fluctuations of the off-diagonal Hamil- 
tonian matrix elements, '~fnmnm = 7i- As stated in the previous section, only fluctuations 
of the strongest interactions, shown in Fig. [H are included. Note that, since 'jnmnm is the 
expectation value of a squared quantity, it is necessarily positive. In the original works 
by Haken, Strobl and Reineker,— interaction fluctuations were already considered, and 
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also a recent paper by Chen and Silbeysi showed that interaction fluctuations can have a 
pronounced effect on exciton dynamics in the dimer. This makes it all the more surprising 
that this effect is commonly neglected in recent studies concerning the FMO complex. The 
time evolution of the population of exciton state s = 3 and the (real part of the) coherence 
between exciton states s = 3 and s = 7 is shown in Fig. [21 clearly showing that their time 
evolution strongly depends on the presence of interaction fluctuations. 




0.1 0.2 0.3 0.4 0.5 



t(ps) 



FIG. 2: Time evolution of the exciton population /J33 (solid lines) and the real part of the exciton 
coherence /O37 (dashed lines), with initial phase 9 = 0, for 70 = 94 cm~^ and 71 = 0, 1,2,5, 10,20 
cm~^. Note the appreciable increase in both transfer rate and decoherence rate with 71. 

Analysis of the temporal evolution of the exciton populations and coherences, in particular 
the ones shown in Fig. |2]and estimating the rates by assuming simple exponential decay and 
looking at times where the oscillations have died out, reveals that the population transfer rate 
and decoherence rate both increase approximately linearly with the interaction fluctuations 
7i. This is to be expected from the general form of the time evolution equations, Eq. |6l 
While these full equations contain a very large amount of terms, all couplings between the 
density matrix elements are proportional to either 70 or 71. The population transfer rate, 
in particular, contains only a small contribution from 70 terms for these parameters, and 
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therefore the population transfer is increased dramatically upon introduction of interaction 
fluctuations. Note that this behavior is consistent with previous results for the dimer,— 
where the exciton population transfer rate for a strongly detuned dimer is to a very large 
extent dominated by the interaction fluctuations. This is analogous to the situation in the 
FMO complex, where the differences between the BChl transition energies are also typically 
large compared to the interactions. More specifically, for a dimer detuned by an energy 
difference Ai? and having an interaction strength J, the mixing angle 9' defined by tan 9' = 
2J/AE is small. It can analytically be shown^i^ that the population transfer rate F between 
the two exciton states is given by F = 70 sin^ 9' + 271 cos^ 9'. However, since the first term is 
typically very small for strongly detuned dimers (i.e., AE ^ J), the population transfer rate 
is fully dominated by the second term as soon as interaction fluctuations are present. An 
analysis of the population evolution shown in Fig. [2] confirms that the population transfer 
rates in FMO are indeed also proportional to 71. Thus, interaction fluctuations lead to 
dramatic increases in population transfer rates, in turn enhancing the efficiency of the FMO 
complex. 

To further illustrate the increase in population transfer rates with increasing interaction 
fluctuations. Fig. [3] shows the harvesting time as a function of the interaction fluctuation 
amplitude 71. In Fig. |3l we consider the energy fluctuation magnitude 70 that optimizes 
the harvesting time, 70 = 94 cm~^, and an increase respectively decrease of one order of 
magnitude. Note that increasing 71 leads to more efficient light harvesting due to a corre- 
sponding increase in population transfer within the FMO complex; this holds for different 
values of the energy fluctuations 70 and for different initial conditions. Additionally, while 
the optimal dephasing value 70 does not change with 71, the optimum becomes considerably 
less deep, making the FMO complex more robust to variations in the dephasing rate 70. For 
very large, in fact unphysically large, interaction fluctuations, the efficiency of the light har- 
vesting complex becomes independent of the energy fluctuation magnitude. This situation 
corresponds to large population transfer rates, dominated by the interaction fluctuations 71, 
and subsequent rapid redistribution of population over all BChl's. The harvesting time will 
in that case converge to r = 7 ps, corresponding to a trapping rate kt = 1 ps^^ from BChl 
3 where at all times l/7th of the total untrapped population will reside. 
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FIG. 3: Dependence of the harvesting time on the interaction fluctuation magnitude 71, for various 
values of the energy fluctuations 79. Notice the decrease in harvesting time with 71 as a result of 
increasing population transfer rates. 

B. Effect of correlated energy-interaction fluctuations 

1. Parameter set 1 

Here, we consider 70 = '^nnnn = 94 cm"^ which corresponds to parameters close to the 



optimal energy transport conditions reported in Ref. llSl The magnitude of the allowed 
correlations is limited by Eqs. |9l to allow for appreciable amounts of correlation, we take 
a large interaction fluctuation value 71 = '^nmnm = 40 cm~^ {n 7^ m), with the same 
provisions as before. In Sec. IIIIB 2\ we discuss the effect of correlated energy and interaction 
fluctuations for a smaller, more realistic value of the interaction fluctuation magnitude 71. 
As stated in Sec. IIIIl we only include correlations that involve the strongest interactions (see 
Fig. [1]) and we focus on correlations between interaction fluctuations and energy fluctuations 
of the BChl's involved in that particular interaction; ■jn'n'nm oc {Sn'n + ^n'm), n ^ m. 

First of all, taking all energy-interaction fluctuation correlations equal in sign and mag- 
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nitude leads to negligible effects. This is consistent with a previous study by Chen and 
Silbey, where it was shown that in a dimer the relevant quantity is the difference between 
the two transition energy-interaction ffuctuation correlations.— While the full time evolu- 
tion equations for the more complicated FMO system are more involved, a similar result 
is observed to hold. This is not surprising, as fluctuation correlations of equal sign would 
imply that a change of interaction would mean that both energies involved change in the 
same direction, leaving the energy difference unchanged. Therefore, we choose the various 
energy-interaction correlations of equal magnitude, but with alternating sign; that is, we 
take 71112 = —72212 = 72223 = —73323 = 72 and so on. There is an ambiguity near BChl 4, 
as it strongly interacts with three other BChl's and it is therefore necessary to choose two 
out of 74434, 74454 and 74474 with equal sign. We take 74434 = 74474 = —74454 = —72; other 
choices result only in inconsequential numerical changes. 

A few observations can be made regarding these plots. First of all, energy-interaction 
fluctuation correlations affect both the energy transfer rates and the coherence decay rates. 
Depending on the sign of the correlations, both increases and decreases of the various rates 
can occur. This can be observed in Fig. HJ where the population of exciton state s = 3 
shows appreciable differences in transfer rate to other exciton states for different correlation 
magnitudes. While this particular exciton state and initial condition shows an increase in 
the transfer rate upon increasing correlations, this is not general; a decrease in transfer 
rate may be obtained for different exciton states and initial conditions. Secondly, it is clear 
that these parameters lead to a quick redistribution of population over all chromophores, 
only showing oscillatory features on a short time scale. Therefore, in the next subsection, 
we consider the same quantities for different parameters, where the effect on population 
oscillations can be observed more clearly. 

It is also possible to investigate the dependence of the harvesting efficiency, here quantified 
through the average trapping time r, on the magnitude of the energy-interaction fluctuation 
correlation 72. This is shown in Fig. [5l the effects are rather limited, with changes of at 
most a few percent. Both the overall effect of accounting for these correlations and the 
dependence on the initial condition, defined by the mixing angle 6', is rather weak. This is 
related to the fact that we are close to the conditions that constitute optimal harvesting 
behavior (i.e., 70 = 94 cm~^), which is corroborated by the harvesting times which are only 
slightly above r = 7 ps. In other words, the bottleneck in the overall harvesting process is 
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FIG. 4: Time evolution of the exciton population P33, 9 = 0, with 70 = 94 cm and 71 = 40 cm . 
The considered correlation values are 72 = —40, —20, 0, 20, 40 cm~^. 

not population transfer between the BChl's, but the trapping of excitations from BChl 3 
into the reaction center. Therefore, the effect of limited changes in the population transfer 
rates have relatively little effect on the overall efficiency. Note that negative values of 72 
may give a small enhancement of the harvesting time; depending on the initial conditions, 
a nonzero optimal choice for the correlations may be made such that light harvesting is at 
its fastest. 

Again, we note that the fluctuations considered here lead to a rapid decay of the co- 
herences. In the next section, we will consider parameters that lead to slower decoherence 
and longer lived population oscillations, in order to evaluate the effect of energy-interaction 
fluctuation correlations on population oscillations. 
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FIG. 5: Dependence of the harvesting time on the magnitude of the energy-interaction fluctuation 
correlations, for various initial conditions and with 70 = 94 cm~^ and 71 = 40 cm""*^. Note that 
negative values of 72 are required for enhancement of the efficiency, and that the optimal correlation 
value depends on the initial condition. 

2. Parameter set 2 

In order to clearly observe the effect of additional energy-interaction fluctuation correla- 
tions on population oscillations, it is necessary to switch to parameters where the decoherence 
occurs on a slower timescale. Here, we consider 70 = 'jnnnn = 10 cm"-*^, and 71 = jnmnm = 2 
cm~^ [n 7^ m) for the interactions shown in Fig. [TJ The same observations as before hold 
with regards to the signs of the correlations: alternating signs are required for observable 
effects and are used throughout this section, while a negative 72 can produce enhancement 
of the transport efficiency. Again, as an example we consider the time evolution of the 
population of the exciton states s = 3 and initial phase ^ = for various values of 72, shown 
in Fig. [61 

It can now be clearly seen that the populations oscillate in time, an effect which is induced 
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FIG. 6: Time evolution of the exciton population p-^^, 9 = and with 70 = 10 cm and 71 = 2 
cm^^. We consider the values 72 = —4, —2, 0, 2, 4 cm^-*^. 

by the coherence that is present between the two exciton states. Indeed, the period of the 
oscillations corresponds to the energy difference between the exciton states, and is identical 
to the coherence oscillation period, which corresponds to Tosc = 130 fs. A clear augmentation 
of the amplitude of the oscillations can be observed when one increases the magnitude of 
72. The introduction of energy-interaction fluctuation correlations can thus enhance the 
population oscillations, which are driven by the presence of coherence. As before, we also 
observe a change in population transfer rates, with in particular a strong suppression of 
population transfer with increasing 72. 

The dependence of the harvesting efficiency on the correlation magnitudes can be calcu- 
lated exactly as before, and is shown in Fig. [71 There is now a stronger dependence on the 
magnitude of the correlations, which is caused by the fact that this parameter set constitutes 
suboptimal conditions, so that in this case changes in the population transfer rates have a 
larger effect on the overall harvesting time. Also, the behavior is almost independent of the 
initial conditions, even more so than for the parameters in the previous section. Negative 
values for 72 are required for enhancement, where the transport is now optimized by choos- 
ing 72 as negative as possible. Positive values of 72 lead to a strong increase in harvesting 
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time; this corresponds to the strong decrease in the relevant population transfer rates, as 
observed in for example Fig. |6l 




FIG. 7: Dependence of the harvesting time on the magnitude of the energy-interaction fluctuation 
correlations, for various initial conditions and with 70 = 10 cm~^ and 71 = 2 cm~^. Note that 
negative values of 72 are required for enhancement of the efficiency. 



C. Effect of correlated interaction-interaction fluctuations 

It is straightforward to include correlations between the various interaction fluctuations. 
As in the previous section, we anticipate that such correlations are strongest between in- 
teractions that have one molecule in common, and in addition we focus on the strongest 
interactions. First of all, it is important to note that while such correlations may very well 
be present, also these are naturally limited in magnitude in order to fulfill the consistency 
conditions Eqs. [9l In particular, we have 

{{Vab - Vacf) OC 7a6afe + lacac " '^.^abac > 0, (16) 
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FIG. 8: Time evolution of the exciton population p^^, = 0, 70 = 94 cm , 71 = 40 cm , and 
72 = 40 cm~^. The considered correlation values are 73 = —30, —15, 0, 15, 30 cm~^. 

where a ^ b,a ^ c. From this, it follows that 

(17) 

Nevertheless, it is instructive to quantify the effects of these additional correlations. The 
conclusions are typically similar to those in the previous section, although the effects are 
typically less pronounced. For brevity, we will not show all the same figures as before, but 
we show a typical result in Fig. |8l Here, we again consider the parameters previously used 
in Sec. HUB 1\ with in addition an energy-interaction fluctuation correlation magnitude 
72 = 40 cm~^ with the sign conventions defined there. First of all, here too alternating signs 
are required to amplify the effect of interaction-interaction fluctuation correlations. We take 
the interaction-interaction fluctuation correlations of the form 71223 = —72334 = 73445 = 73 
et cetera. Again there is some ambiguity in the deflnition of the signs for the correlations 
around BChl 4, however, since the details of those choices hardly matter for the initial 
conditions we have chosen, we will not go in to this in more detail. We show time evolution 
plots only for 6* = 0, behavior for other parameters is not significantly different. First of 
all, it is clear that interaction-interaction fluctuation correlations, too, lead to changes in 
the population transfer and dephasing rates, as exemplifled by the population evolution of 
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FIG. 9: Dependence of the harvesting time on the magnitude 73 of the interaction-interaction 
fluctuation correlations, for various initial conditions. We have chosen 70 = 94 cm^"*^, 71 = 40 
cm~^ and 72=40 cm~^. The different offsets of the curves are consistent with the observed 6- 
dependence of the harvesting time in Fig. [5l 

exciton state 3 shown in Fig. |8l In addition, while there is a small change in the size of 
the populations, the amplitude of the oscillations is influenced by the interaction-interaction 
fluctuation correlation magnitude to only a small extent. These conclusions are observed to 
hold for different sets of parameters and initial conditions. 

Fig. |9] shows the dependence of the harvesting time on the magnitude of the interaction- 
interaction fluctuation correlations. Again, there is an optimum for some nonzero value of 
the correlation magnitude 73; the exact value of the optimal magnitude for the correlations 
depends on the initial conditions and on the choice for the energy-interaction correlations. 
The behavior for other values of the initial phase factor 9 is very similar, implying that 
the transfer rates and other relevant parameters for trapping will only depend weakly on 
possible interaction-interaction fluctuation correlations. While Fig. |9] shows that the effect 
of interaction-interaction fluctuation correlations is roughly of the same magnitude as the 
effect of transition energy-interaction fluctuation correlations, this is only the case if the 
former occurs in conjunction with the latter. The changes induced by interaction- interaction 
fluctuation correlations for the case of 72 = cm~^ are considerably smaller. Finally, as 
was previously the case in Sec. IIIIBl a more suboptimal choice of parameters leads to a 
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larger effect of changes in transfer rates, and will thus also amplify the effect of including 
interaction-interaction fluctuation correlations to an extent. 

IV. CONCLUSIONS 

Since in principle environmental changes will induce fluctuations of both transition en- 
ergies and interactions that should be correlated to some extent, we have performed a 
study of the effect of such fluctuations and their possible correlations on the energy transfer 
and excitation dynamics in photosynthetic complexes, specifically focusing on the Fenna- 
Matthews- Olson complex as a model system. The inclusion of interaction fluctuations and 
correlations between transition energy fluctuations and interaction fluctuations, and between 
different interaction fluctuations, has been studied by a straightforward generalization of the 
Haken-Strobl-Reineker model. These additional correlations have been shown to be natu- 
rally limited in magnitude, in order to correspond to a physically consistent system. 

Interaction fluctuations on the excitation dynamics in photosynthetic complexes can in 
general not be neglected. Not only will these in general occur, but it has been shown that 
even small values of interaction fluctuations can already lead to an appreciable increase in 
transfer rates and a corresponding increase in efficiency of the FMO complex. In particu- 
lar, for suboptimal values of the energy fluctuations, interaction fluctuations can lead to a 
considerable optimization of the light harvesting process, even to such an extent that the 
transfer rates and the eventual efficiency is dominated by interaction fluctuations and not by 
energy fluctuations. The presence of interaction fluctuations may thus make the efficiency 
of the FMO complex more robust to changes in the dephasing rate. Generally, for various 
initial conditions and energy fluctuation magnitudes, an increase in the efficiency of the 
FMO complex is observed with increasing interaction fluctuation amplitude. 

When including correlated transition energy fluctuations and interaction fluctuations, it 
is first of all possible to observe somewhat enhanced population oscillations, driven by the 
coherence between the initially excited exciton states. Depending on the signs of the various 
correlations, a suppression of the oscillations can also occur. In addition, one sees that 
the transfer rates between exciton states is modified, which in turn leads to accompanying 
changes in the overall efficiency of the photosynthetic complex. The overall effects on the 
efficiency are limited in scope when one considers parameter values that approximately 
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correspond to an optimal functioning of the light harvesting complex. In that case, not 
the excitation transfer but the trapping into the reaction center is the bottleneck, so that 
changes in transfer rates will only change the harvesting time by up to a few percent. For less 
optimal conditions, the effects can be considerably larger, as the excitation transfer processes 
play a larger role in determining the overall harvesting time. Typically, a small increase in 
the light harvesting efficiency can be obtained by using a nonzero amount of correlation 
between the aforementioned fluctuations. By the same approach, one can quantify how 
correlations between the various interaction fluctuations modify previously obtained results. 
The effect of these additional correlations turns out to be qualitatively very similar to that of 
transition energy-interaction fluctuation correlations, and are quantitatively of a comparable 
magnitude. The net effect of such correlations is limited to small changes in the scattering 
rate and resultant changes in the overall efficiency. 

This study thus suggests that one should not neglect interaction fluctuations in models 
describing the energy transfer in photosynthetic complexes. Correlations may also play a 
relevant role in the excitation dynamics, and this may (but does not necessarily) translate 
into appreciable changes in the efficiency. 
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Appendix A: Exciton states and energies in the FMO complex 

The exciton states are found by diagonalization of the Hamiltonian matrix (fT2|) : their 
coefficients and energies are given in the table below. 
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